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Abstract 

Background: Hypermethylation of CpG islands in tumor suppressor gene plays an important role in carcinogenesis. Many 
studies have demonstrated that hypermethylation in promoter region of RARfi gene could be found with high prevalence in 
tumor tissue and autologous controls such as corresponding non-tumor lung tissue, sputum and plasma of the NSCLC 
patients. But with the small number subjects included in the individual studie, the statistical power is limited. Accordingly, 
we performed this meta-analysis to further asses the relationship of methylation prevalence between the cancer tissue and 
atuologous controls (corresponding non-tumor lung tissue, sputum and plasma). 

Methods: The published articles about RARp gene promoter hypermethyltion were identified using a systematic search 
strategy in PubMed, EMBASE and CNKI databases. The pooled odds ratio (OR) of RARP promoter methylation in lung cancer 
tissue versus autologous controls were calculated. 

Results: Finally, eleven articles, including 1347 tumor tissue samples and 1137 autologous controls were included in this 
meta-analysis. The pooled odds ratio of RARp promoter methylation in cancer tissue was 3.60 (95%CI: 2.46-5.27) compared 
to autologous controls with random-effect model. Strong and significant correlation between tumor tissue and autologous 
controls of RARP gene promoter hypermethylation prevalence across studies (Correlation coefficient 0.53) was found. 

Conclusion: RARP promoter methylation may play an important role in carcinogenesis of the NSCLC. With significant 
methylation prevalence correlation between tumor tissue and autologous of this gene, methylation detection may be a 
potential method for searching biomarker for NSCLC. 

Citation: Hua F, Fang N, Li X, Zhu S, Zhang W, et al. (2014) A Meta-Analysis of the Relationship Between RARft Gene Promoter Methylation and Non-Small Cell 
Lung Cancer. PLoS ONE 9(5): e96163. doi:10.1371/journal.pone.0096163 

Editor: Qing-Yi Wei, Duke Cancer Institute, United States of America 

Received November 27, 2013; Accepted April 4, 2014; Published May 5, 2014 

Copyright: © 201 4 Hua et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This project was supported by the grants from National Natural Science Foundation of China (to XBL) (No. 81302002). Science and technology 
foundation of Tianjin Municipal Health Bureau (to JDG) (2010KZ040). The funders had no role in study design, data collection and analysis, decision to publish, or 
preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: jundonggu@aliyun.com (J-DG); Weisan_zhang@yeah.net (W-SZ) 
9 These authors contributed equally to this work. 



Introduction 

Lung cancer, accounting for 100 million deaths and 120 million 
incidence per year world-wide, is the leading cause of cancer 
related death[l]. Generally, lung cancer is divided in to non-small 
cell lung cancer (NSCLC) and small-cell lung cancer according to 
the pathology. The prognosis of NSCLC, accounting for 80-85% 
of the lung cancer, is poor for lack of effective early diagnosis 
methods [2]. Recently, CpG island methylation of tumor suppres- 
sor genes promoters have been found to be important for 
carcinogenesis in many kinds of carcinomas including NSCLC [3]. 
The DNA methylation is a biochemical process involving the 
addition of a methyl group to thecytosine or adenine DNA 
nucleotides. DNA methylation may affect the transcription of the 
tumor suppressor genes in two ways. Firstly, the methylated DNA 
may physically impede the binding of transcriptional proteins to 
the gene itsel [4] . Secondly, may be with more important, 
methylated DNA may be bound by proteins known as methyl- 



CpG-binding domain proteins (MBDs)[5]. MBD proteins then 
recruit additional proteins to the locus, such as histone 
deacetylases and other chromatin remodeling proteins that can 
modify histones, thereby forming compact, inactive chromatin, 
termed heterochromatin [6] . 

RARP gene is a typical tumor suppressor gene encoding retinoic 
acid receptor beta, a member of the thyroid-steroid hormone 
receptor subfamily, which belongs to nuclear receptor superfamily. 
The protein encoded by this gene could bind retinoic acid, the 
biologically active form of vitamin A which mediates cellular 
signallings in embryonic morphogenesis, cell growth and differen- 
tiation [7]. Recent studies demonstrated that hypermethylation in 
promoter region of RARfi gene could be found with high 
prevalence in tumor tissue and autologous controls such as 
corresponding non-tumor lung tissues (CNTLT), sputum and 
plasma of the NSCLC patients. But with the small number 
subjects included in the individual studies, the statistical power is 
limited. Thus, a meta-analysis about the relationship of RARP 
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promoter methylation between lung cancer tissue and autologous 
controls was done in order to better identify the correlation of 
methylation status between cancer tissue and autologous samples. 

Methods 

Studies identification 

The selection procedure of the relevant articles was demon- 
strated in figure 1 . Papers about RARfl gene promoter CpG island 
hypermethylation in NSCLC, published before February 2013, 
were identified by searching Pubmed, EMBASE and CNKI 
databases. The "Non-small cell lung carcinoma/NSCLC" 
"methylaiton/hypermethlation" and " RAR/i / 'RAR-beta" were used 
as the MeSH and free text word when searching the databases. 
The inclusion criteria: non-small cell lung cancer patients with 
histology or cytology conformation; The methylation was detected 
by methylation-specifec PCR (MSP) or quantitative MSP (q-MSP); 
The methylation rate in the cancer tissue and autologous controls 
would be drawn in the original studies included in this meta- 
analysis. All potentially relevant articles were assessed in full-text 
paper and all references of included articles were scanned for 
additional analysis according to the Cochrane handbook for 
systematic review [8]. 



Data extraction 

Information about the general characteristics and methylation 
prevalence of the original included studies were retrieved. For 
general characteristics: The first author, year of publication, 
journals, race of the included subjects, mean or median age of the 
included patients, methylation detection methods were all 
recorded. For methylation prevalence in cancer tissue and 
autologous controls: the RARp promoter methylation rate or 
number of methylated (M) and unmethylated (U) samples in 
cancer tissue and autologous controls were all retrieved from the 
included studies. The methylation prevalence (MP) was calculated 
according to the formula PM = M/(M+U) x 100%. All of the 
information details were extracted by two reviewers (FH and NZF) 
and then checked by the third reviewer (JDG). 

Statistical analysis 

The odds ratio (OR) for each individual study was calculated by 
the formula OR = (a/b)/ (c/ d). (a = number of methylated samples 
in tumor tissue b = number of unmethylated samples in tumor 
tissue; c = number of methylated in control tissue and d = number 
of unmethylated in control tissue.) The OR and its 95% 
confidence intervals (CI) of RARfl promoter hypermethylation in 
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Figure 1. Flowchart of the literature search strategy (Data from some studies was used more than once). 

doi:10.1371/journal.pone.0096163.g001 
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Figure 2. Forest plot of pooled OR for RAR/Sgene promoter methylation in cancer tissue versus autologous controls. 

doi:1 0.1 371 /journal.pone.00961 63.g002 



cancer tissue compared to autologous controls (CNTLT, sputum 
and plasma) for the included studies was pooled by STATA/SE 
1 1.0(StataCorp LP, http://www.stata.com)statistic software. Sta- 
tistical heterogeneity among articles included in this meta-analysis 
was evaluated by I 2 [9] .Without significant heterogeneity (I 2 > 
50%), the fixed-effect method was used to pool the data. Inversely, 
with significant heterogeneity among studies, random-effect 
method (Dersimonian-Laird method) was taken to pool the OR. 
The methylation prevalence correlation between cancer tissue and 
autologous clinical control (CNTLT, sputum and plasma) was 
evaluated by Spearman's rank correlation test. And the publica- 
tion bias was assessed by Begg's funnel plot and Egger's test [10]. 

Results 

General characteristics of included studies 

Finally, eleven articles [1 1—21], including 1347 tumor tissue 
samples and 1137 autologous controls, were included in this meta- 
analysis. (Figure 1). Of the 11 articles included in this meta- 
analysis, 6 were conducted in Chinese mainland, 2 were in 
Taiwan, 2 in USA and 1 in Italy. The median methylation 
prevalence of RARji gene promoter was 47.56% and 17.50% in 
cancer tissue and autologous controls, respectively (table SI). The 
general characteristics of the included articles were showed in 
table 1. 



Pooled results of this meta-analysis 

Significant heterogeneity was found across the included studies 
(I 2 = 72.8%). Accordingly, the random-effect model was taken in 
pooling the OR. In this meta-analysis, the pooled odds ratio of 
RARfi promoter methylation in cancer tissue was 3.60 (95%CI: 
2.46-5.27, p = 0.000) versus autologous controls with random- 
effect model (Figure 2). The sensitivity analysis was also done by 
omitting a single study under the random-effect model. But no 
significant changes of OR and it 95%CI was found in the 
sensitivity analysis which demonstrated that the pooled odds ratio 
was not sensitive to a single study. 

Subgroup analysis 

With significant heterogeneity across the included studies 
(I 2 = 72.8%), the subgroup analysis of this meta-analysis was 
performed according to race (East-asia, Caucasian) and control 
types (corresponding non-tumor lung tissues, sputum and plasma). 
In the subgroup analysis, the significant odds of the RARP 
promoter methylation in tumor tissue was only changed when 
comparing to sputum (OR= 1.33, 95%CI: 0.98-1.80, P=0.120) 
and in subjects with Caucasian ethnicity (OR = 2.38, 95%CI: 
0.71-7.92, P= 0.159, Table 2). However, the conclusion of the 
sputum control and Caucasian ethnicity subgroup should be 
interpreted with caution as only a two and one articles with small 
subjects were included in these subgroups analysis. 
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Figure 3. The distribution of methylation rate for each included 
study (the paired tissue were in the same colour). 

doi:10.1371/journal.pone.0096163.g003 

Methylation prevalence correlation 

The methylation prevalence was much higher in tumor tissue 
(meidan methylation prevalence 47.56%) compared to autologous 
controls (meidan methylation prevalence 17.50%) of RARji gene 
(Figure 3). And we also found significant RARfS gene promoter 
hypermethylation prevalence correlation between tumor tissue 
and autologous controls among the included studies (Correlation 
coefficient 0.53) (Figure 4). It means that the more methylation 
percentage in tumor tissue and the more methylation percentage 
in autologous controls. 

Publication bias 

No statistical significant publication bias was detected by 
Begger's funnel plot and Egger's tests (t=0.70, p = 0.500) 
(Figure 5). 

Discussion 

Lung cancer is leading cause of cancer related mortality world- 
wide with 100 million deaths and 120 million diagnoses in year 
2009 [22]. Without effective early diagnosis method, the prognosis 
of NSCLC is relative poor with a 5-year survival rate of less than 
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Figure 4. Correlation of RAR// gene promoter methylation 
between tumor tissue and autolougs controls. 

doi:10.1371/journal.pone.0096163.g004 
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Figure 5. Begg's funnel plot for assessment of publication bias. 

doi:10.1371/journal.pone.0096163.g005 



1 5 % [2] . CpG island methylation of tumour suppressor genes 
promoters have been found to be important for carcinogenesis and 
early events in many kinds of carcinomas including NSCLC. 
RARfi gene is a typical anti-tumor gene encoding retinoic acid 
receptor beta, a member of the thyroid-steroid hormone receptor 
subfamily of nuclear receptor superfamily. The protein encoded by 
this gene could bind retinoic acid, the biologically active form of 
vitamin A which mediates cellular signallings in embryonic 
morphogenesis, cell growth and differentiation. This gene was 
always been silence by its promoter DNA hypermethylation on its 
promoter region which can physically impede the binding of 
transcriptional proteins to the gene itself. 

The relationship between cancer tissue and autologous controls 
such as CNTLT, plasma and sputum for this gene was not clear 
with small subjects in the individual studies. So, we performed this 
meta-analysis, which can strength the statistical power to quantify 
the hypermethylation-disease association by pooling data from 
open published articles. In this meta-analysis, we finally included 
eleven articles with 1347 tumor tissue samples and 1137 
autologous controls to pool the odds ratio of methylation 
prevelance. The pooled odds ratio of RARfi promoter methylation 
in cancer tissue was 3.60 (95%CI: 2.46-5.27) compared to 
autologous controls with random-effect model. The subgroup 
analysis of this meta-analysis was performed according to race 
(East-asian, Caucasian) and control types. In the subgroup 
analysis, the significant odds of the RARfi promoter methylation 
in tumor tissue was only changed when comparing to sputum 
(p = 0.120) and in subjects with Caucasian ethnicity (p = 0.159). 
However, the conclusion of the sputum control and Caucasian 
ethnicity subgroup should be interpreted with caution as only a 
two and one articles with small subjects were included in these 
subgroup analyses. We also found strong and significant correla- 



tion between tumor tissue and autologous controls of RARfi gene 
promoter hypermethylation prevalence across studies (r = 0.53). 
The correlation showed that, the higher methylation prevalence in 
plasma/sputum samples, the higher methylation prevalence could 
be found in cancer tissue in patients with NSCLC. And this 
indicated that detection methylation status in plasma or sputum 
could be a useful assay for diagnosis of NSCLC without or with 
mirror invasion. Nevertheless, in our meta-analysis, all of samples 
(tumor tissue, non-tumor lung tissue sputum and plasma) are all 
come from the patients with confirmed lung cancer. There was no 
relative healthy person controls. So, the lung cancer screening 
sensitivity and specificity was not possible to calculate according to 
the Bayes theorem. 

Some limitations of this meta-analysis are need for further 
consideration. First, significant heterogeneity was found in this 
meta-analysis (I 2 = 72.8%), which could decrease the statistical 
power of this study. Second, the co-variate methylaiton status of 
different tumor suppressor genes may also be linked and interact 
with each other, suggesting methylation analysis of a single gene 
may be not enough [2 3]. 

In conclusion, this meta-analysis showed RARfi promoter 
methylation was much higher in tumor tissue compared to 
autologous controls, which indicates promoter hypermethylation 
of tumor suppressor gene may play an important role in 
carcinogenesis of the NSCLC. With significant methylation 
prevalence correlation, methylation detection may be a potential 
stragegy for searching biomarker for NSCLC [24] . 

Supporting Information 

Table SI The methylation percentage in each included study 
(DOC) 
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